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4. Problemas resueltos 


4.1. Ecuaciones diferenciales ordinarias de primer orden 


Ejercicio 4.1. Resolver las siguientes ecuaciones diferenciales ordinarias de primer orden. 

A. Variables separadas 

b) (1 +e x )yy' = e x . Hallar la solución que pasa por (0,1). 

B. Homogéneas 

C. Lineales 

h) xy’ + 4y = x 3 — x. 

D. Diferenciales exactas 

e) (x 3 +xy 2 ) dx+{x 2 y+y 3 ) dy = 0. 


Resolución. 

A. Variables separadas 

b) En primer lugar buscamos la solución general de la ecuación diferencial. Separando las variables e 
integrando, 


6 X f f 6 X 

(1 + e*)y dy = e x dx => y dy = ——- dx => ydy= — —- dx 

í+e x J J 1 + e x 

y 2 

— = ln(l +e x ) +c, 

de donde obtenemos 

y 2 = 21n(l +e x ) +C => y = ± </21n(l +e x ) +C. 

Para obtener la solución que particular que pasa por (0,1) consideramos la solución positiva de la ecuación 
diferencial, esto es, y — \J2 1n(1 +e x ) +C. Como ésta ha de pasar por (0,1), se debe tener y(0) = 1. Por tanto: 

y(0) = \/21n(l +e°)+C = 1 v / 2lñ2 + C=l ^ 21n2 + C = l 

C = 1 — 21n2. 
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Luego, la solución particular buscada viene dada por: 


y = -J 21n(l+e JC ) + l-21n2 


y = ^/2[ln(l +e x ) — ln 2] + 1 


y = 


1 +21n 


1 +e x 
2 


B. Homogéneas 

d) Nótese que la ecuación diferencial puede ser escrita como 


dy = y 

dx x 


Haciendo el cambio de variable z = - ó, lo que es lo mismo, y 

x 


= xz, se obtiene 


dy dz 

dx " dx 


Encontramos así que 


dz . 

z+x — = 1 —z 
dx 


dz i o 

x — = 1 — 2z 
dx 

1 (-2) , dx 
a , dz — 

2 1 — 2z x 

ln(l—2 z) _ 2 = lnx + lnC 


dz dx 
1—2 z x 

— ^ ln( 1 — 2z] = lnx + c 
ln(l — 2z) - z = InCx 


=> (1—2 z) 1 = Cx 


(1 — 2z) = Cx“ 2 


2z = 1 - ^ 


y 

X 


C 

y-2 


Z = 2 


y = 


c 

X 2 

C 


Calculamos ahora la solución particular que verifica y( 1) = 1: 


y(l)= 1 --C=\ 



1 

2 ' 


Por tanto, la solución particular buscada viene dada por 


x 


y = 


i 

2x 


C. Lineales 
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h) Dividiendo por x ambos miembros de la ecuación diferencial resulta 


/ . 4 2 i 

y + - y = x — 1. 
x 


Multiplicamos la ecuación por el factor integrante e-f t dx = e 41nx = e ln ' 4 =x 4 , con lo que obtenemos 


x 4 y' + 4x 3 y = x 6 —x 4 => — ( x 4 y ) = x 6 — 


.6 v 4 


dx 


X' X 5 


y = /(x — x 4 )dx 


_r x C 


A=tt-^+c => y= —-- + ^7. 


7 5 


7 5 x 4 ' 


D. Diferenciales exactas 


e) Escribiendo M(x,y ) = x 3 + xy _ y N(x,y) = x 2 y +y , se cumple que 


dM d(x 3 +xy 2 ) . <9(x 2 y + y 3 ) <9A, 

=- -5 -= 2 xy =-^-= - 3 —(x,y). 

ay ay dx dx 


Al coincidir ambas expresiones se trata, en efecto, de una ecuación exacta. Por tanto. 


4 2 2 

i(x,y) = j M ( x ,y) dx + c(y) = J (x 3 +xy 2 ) dx + C(y) = + -2- +C(y). 


Para calcular C(y), 


dy 


(x,y)=x 1 y+C'(y)=N(x,y)=x 2 y + y 3 => C(y)=y 3 


C(y) = ^+C. 


Consecuentemente, la solución general será 


.4 .2 ,2 ,4 

w(x,y)=C => ^ ^ ^ = C (C constante), 


que podemos reescribir como 


- 4 + 2x 2 y 2 + y 4 = C (C constante), 


o bien 


(x 2 +y 2 ) 2 = C (C constante). 


□ 
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4.2. Aplicaciones: problemas de desintegración radiactiva 


Ejercicio 4.2. Inicialmente había 100 miligramos de una sustancia radiactiva. Después de 6 horas su masa 
disminuyó en un 3 %. Si en un instante cualquiera la rapidez de desintegración es proporcional a la cantidad 
de sustancia presente, determinar la cantidad que queda después de 24 horas. 


RESOLUCIÓN. Inicialmente tenemos 100 mg de sustancia radiactiva. Si C(t) denota la cantidad de sustancia 
radiactiva en el instante t , sabemos que al cabo de t = 6 h quedan 

C( 6 ) = 100-3 = 97 gr 

de esta sustancia. La rapidez de desintegración es proporcional a la cantidad de sustancia presente, esto es: 

f-. 

siendo k la constante de proporcionalidad. Como vimos en el desarrollo teórico, tal ecuación admite por solu¬ 
ción C{t) = A e kl . donde Ay k son constantes a determinar. Puesto que en el instante inicial t = 0 contamos 
con 100 mg de sustancia, 

C(0) = A e° = 100 => A = 100. 

En el instante t = 6 quedan 97 gr; luego, 

97 97 1 97 

C( 6 ) = 100= 97 =► e<*=- =► 6 * = ln- => 

En conclusión, la cantidad de sustancia radiactiva en el instante t es 

C(t) = lOOe® 1 "™. 

Por tanto, la cantidad remanente transcurridas 24 h es 

C(24) = lOOe^ 111 ™ =100e 41n ™ = 100 e ~ 012 ~ 88.5 mg. 


□ 
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4.3. Aplicaciones: problemas de mezclas 

Ejercicio 4.3. Un tanque está lleno de 100 litros de agua en los que se ha disuelto 20 kilogramos de sal. 
Otra mezcla que contiene 1 kilogramo de sal por litro es bombeada al tanque a razón de 7 litros por minuto. 
La solución mezclada es bombeada hacia el exterior a razón de 8 litros por minuto. Determinar la función 
que da la cantidad de sal en cada instante. ¿Se vaciará totalmente el tanque? 


RESOLUCIÓN. Conforme a la notación utilizada en el desarrollo teórico, el enunciado del problema propor¬ 
ciona los siguientes datos: 


A = 20 kg, a = 1 kg/L, 

Vo = 100 L, vi = 7 L/min, V 2 = 8 L/min. 

Por tanto, la ecuación diferencial que modeliza la cantidad de sal en el tanque en cualquier instante viene dada 
por 


y'it)- 


8 


100+ (7 — 8)f 


y(t) = 7-i, 


esto es, 


yw+ iooU^ = 7 ' 


La ecuación anterior admite como factor integrante 


e í IT = e ~ s f wnh dl — ¿,-81n(100—f) _ ^l+lOO-r) 8 _ _ 


= e 


= e 


(100-f) 8 ' 


Multiplicando ambos miembros de la ecuación por el factor integrante, 


1 „ x 8 . , 7 

■y (0 + /.AA y(t) = 


(100-f) 8 (100-í) 9 (100-f) 8 ’ 


es decir: 


Integrando la expresión anterior. 


d 

Jt [(100-f) 8 


y(t) 


(100-f) 8 ' 


( m - t )* Át) I (ioo-f ) 8 dt ^ 7 \í 


dt — / , „„ - — + C, 


(100-f) 8 (100-f) 7 
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de modo que la solución general de la ecuación diferencial es 


, t \ í 

y(t) = (100-/) +(100-f) 8 C= (100-í)+ 100 8 (^1- — J 


C. 


Para hallar C tenemos en cuenta que la concentración inicial es A = 20: 


A = y(0) = 20 => 20 = 100+ 100 8 C => 100 8 C =-80. 


En conclusión, la cantidad de sal presente en el tanque en cada instante es 


y(t) = (100 — f)—80^1 — — ) • 


100/ 


Para averiguar si el tanque se vaciará totalmente, determinaremos el tiempo en que la concentración se 
anula, esto es: 


y(t) = o =► (ioo-/)-8o(i-^) 8 = o => (loo-/) 


1-80 


(100-í 

100 8 




= 0 . 


La ecuación anterior admite dos soluciones: 


100-í = 0 


t = 100 min; 


(100-í) 7 

^Soo 8 ^ 0 


= 100-f 


100 8 , .7 7/ 100 8 

— = (100-1) * VA» 

7 /100 8 

t = 100- y ~ 100- 103.23 = -3.23 min. 


La solución negativa carece de sentido en el contexto del problema. Por tanto, la concentración es cero para 
t = 100 min, que es cuando se vaciará el tanque. 

Nótese que aunque éste se vacíe siempre seguirá entrando agua salada, de manera que a partir del instante 
t = 100 min la concentración de sal en cada instante será la de la mezcla entrante, a saber, 1 kg/L. □ 


4.4. Aplicaciones: problemas de temperatura 

Ejercicio 4.4. Cuando un objeto absorbe calor del medio que lo rodea sigue la Ley de Newton. Una pequeña 
barra de metal, cuya temperatura inicial es de 20 °C, se deja caer en un recipiente con agua hirviendo. 
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Calcular el tiempo que dicha barra tardará en alcanzar los 90 °C, si se sabe que su temperatura aumentó 2 
°C en un segundo. ¿Cuánto tardará la barra en alcanzar los 98 °C? 


RESOLUCIÓN. La Ley de Newton expresa que la rapidez con que se enfría un objeto es proporcional a 
la diferencia entre su temperatura y la temperatura ambiente. La ecuación diferencial que modeliza dicho 
fenómeno es T'{t) = K\T{t) — T a \ . cuya solución general (veáse el desarrollo teórico) es 

T(t) = T a +Ce Kt . 

La temperatura ambiente en este caso es T a = 100, mientras que la temperatura inicial es T (0) = 20. Por 
tanto, 

L(0) = 100 + Ce*° = 100 + C = 20 => C = —80. 

Como la temperatura aumentó 2 °C en 1 s encontramos que T(l) = 22. Así, 

7 © 

r(l) = 100-80^ = 22 => 78 = 80 e K => K = ln—. 

80 

Por tanto, la temperatura en cualquier instante t es 

T{t) = 100 —80e íln ®. 


Para calcular el tiempo que tarda la barra en alcanzar 90 °C resolvemos la ecuación T (/) = 90: 


100 — 80 e íln sD =90 


10 = 80 e ,ln w 


i=e íln m 
8 


t = 


lll; 


l n 7l 
111 80 


~ 82.1 s. 


Similarmente, para calcular el tiempo que tarda en alcanzar 98 °C resolvemos la ecuación T(t) = 98: 


100 — 80 e íln ® =98 


2 = 80 e íln ® 


40 


= e íln S5 


=> t = 


ln 40 


ln 


~ 145.7 s. 


□ 
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4.5. Aplicaciones: miscelánea 

Ejercicio 4.5. La rapidez con que cierto medicamento se disemina en el flujo sanguíneo se rige por la 
ecuación diferencial 

d -^=A~BX, X(0)=0, 
dt 

donde A y B son constantes positivas. La función X(t) describe la concentración del medicamento en el 
flujo sanguíneo en un instante cualquiera t. Encontrar el valor límite de X cuando t —> °°. ¿Cuánto tarda la 
concentración en alcanzar la mitad de este valor límite? 

Resolución. La ecuación diferencial 


dX 

dt 


= A-BX, X{0)=0, 


es de variable separada (nótese que también es de tipo lineal). Resolvemos por el procedimiento habitual: 


dX 

A-BX 


= dt 


" dX 
A-BX 
1 


= dt 


-ln (A-BX)=t + C 

B 


=> A — BX = e 


— ¿,—Bt+C 


1 f ( -B) 


dX = dt 


Bj A-BX 
ln (A-BX) = -Bt+C 

BX =A — Ke~ Bt 


=> X=--Ke~ Bt . 

B 


Ahora imponemos la condición inicial: 


X(0) =0 




A 

K = —. 
B 


Por tanto, la concentración en el instante t viene dada por 


X(t)=-( l-e-™). 


El valor límite de X ( t ) cuando t —> °° es: 

lím X(t) = — (1 — líme~ Bí ) = —, 
t^OO W B V t^OO ' B 

puesto que la constante B es positiva. Para saber cuánto tiempo se tarda en alcanzar la mitad de ese valor límite. 
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resolvemos 


m 


1 A 

2 B 


A 

B 


( 1 - e- Bt ) 


1 A 

2 B 


1 — e 


Bt 


1 

2 


-Bt 


1 

2 



Bt = ln - 
2 


— In2 


□ 
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